Long-range/short-range separation of the electron-electron interaction in density 

functional theory 
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By splitting the Coulomb interaction into long-range and short-range components, we decompose 
the energy of a quantum electronic system into long-range and short-range contributions. We show 
that the long-range part of the energy can be efficiently calculated by traditional wave function 
methods, while the short-range part can be handled by a density functional. The analysis of this 
functional with respect to the range of the associated interaction reveals that, in the limit of a very 
short-range interaction, the short-range exchange-correlation energy can be expressed as a simple 
local functional of the on-top pair density and its first derivatives. This provides an explanation for 
the accuracy of the local density approximation (LDA) for the short-range functional. Moreover, 
this analysis leads also to new simple approximations for the short-range exchange and correlation 
energies improving the LDA. 



I. INTRODUCTION 

In the Kohn-Sham (KS) approach Q of density func- 
tional theory (DFT) Q of inhomogeneous electronic sys- 
tems, the central quantity is the unknown exchange- 
correlation energy functional E xc [n] which encompasses 
all the many-body effects. The vast majority of approx- 
imations for this functional are based on the original lo- 
cal density approximation (LDA) 2], an approximation 
that turned out to be more accurate and reliable than 
expected and rather difficult to improve in a systematic 
way @- 

Actually, it has been realized for long, with the wave- 
vector analysis of E xc [n] by Langreth and Perdew 
that the LDA describes accurately (but not exactly in 
general 0) short wavelength density fluctuations, but 
is inadequate for long wavelength fluctuations. A dual 
analysis in real space of the exchange-correlation en- 
ergy (see, e.g., Ref. 6) leads to the same conclusion 
that the LDA is accurate at small interelectronic dis- 
tances but fails at large distances. This observation lead 
to the development of the first gradient corrected func- 
tional 0, B M ES IH 111 with the basic objective to 
cure the wrong long-range contribution to the exchange- 
correlation energy of the LDA. 

However, describing accurately the non-local correla- 
tion effects arising from the long-range character of the 
Coulomb interaction by (semi)local density functional 
approximations still seems out of reach. This idea in 
mind, it is has been proposed to use a density functional 
approximation only for the short-range part of the elec- 
tronic energy, and treating the lon g-ra nge part by more 
appropriate many-body methods (Dll Il4 ITEl llri Il7j . 
This approach was somehow inspired from early calcula- 
tions of the correlation energy of the uniform electron gas 



based on a separate treatment of long-range and short- 
range contributions (see, e.g., Refs. ll&ll&l2(]|) . 

For atoms and molecules, this approach leads to a rig- 
orous method for combining traditional ab initio wave 
function calculations with DFT 0, 122^ . and we recall 
now the formalism. The Coulomb electron-electron in- 
teraction is decomposed as 



i=t£(r)+i&(r), 



(1) 



where v£ e (r) is a long-range interaction, vg e (r) is the 
complement short-range interaction and /i is a parameter 
controlling the separation. The universal functional |2^| 
F[n] = min^^„(*|T + V ee \^), where f is the kinetic 
energy operator and V ee = J2i<j ^l r ij lSl the Coulomb 
interaction operator, can then be decomposed as 



F[n] = F»[n] + F*[r 



(2) 



where F^ [n] is the universal functional corresponding to 



the long-range interaction V£, = 



i<3 



F^[n} = min(*|T + V^|*) 



(3) 



and F M [n] = F[n]— F^n] is by definition the complement 
(short-range) part. 

The exact ground-state energy of an electronic system 
in the external local nuclei-electron potential v ne (r) can 
be written using this short-range functional F^n] via 
application of the variational principle 

E = min ^F fl [n]+F^[n\ + J n(r)v ne (r)dr} 

»4»|<t|* + <!s|t> + fM + /«.W*.W*} 



(*"|f + Vi i e \^) + F"[n^}+ / 7i^(r)tw(r)dr 
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(4) 
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In this equation, ^ M is given by the Euler-Lagrange equa- 
tion 



(5) 



where H» = T + V£ + V». Thus, 1 is the ground- 
state multi-determinantal wave function of a partially 
interacting system with interaction V£ and external lo- 
cal potential V 1 * = J2i v>1 ( r i) where v^(r) = v ne (r) + 
8F^[n^]/ 5n{r). By virtue of the Hohenberg-Kohn the- 
orem 0, u M (r) is the unique potential (up to an additive 
constant) which insures that the ground-state density of 
this fictitious system is identical to the ground-state den- 
sity of the physical system = n. 

The short-range functional F M [n] is further decom- 
posed as 



F»[n] = U^r 



EH c [ 



(6) 



where C/ M [«] = 1/2 JJ n(T 1 )n{Y 1 )v^ e {r 12 )dv 1 dv2 is the 
complement short-range Hartree energy functional and 
E£ c [n] is the unknown complement short-range exchange- 
correlation energy functional. Once an approximation is 
chosen for E£ c [n], the wave function can be computed 
by solving self-consistently Eq. J5J using ab initio wave 
function methods like configuration interaction (CI) or 
multi-configurational self-consistency field (MCSCF) and 
the total energy is calculated according to the last line 
of Eq. (0J. Notice also that it is possible to separate the 
functional 25£ c [n] into short-range exchange and correla- 
tion contributions, E!^ c [n} = E!^[n] + E^ [n], with 



and 



EH[n] = ^\VfM-U^[n], 



Eg[n]=P"[n]-($\V£\$) 



(7) 



(8) 



where $ is the KS determinant and V£ = V ee — V£, is 
the complement short-range interaction operator. 

A simple approximation for E ^Jn] is the LDA associ- 
ated to the modified interaction |17ll24j 



been efficiently implemented into the quantum chem- 
istry program Dalton, performing the coupling of MC- 
SCF wave function calculations with DFT [25|. For a 
reasonable long-range/short-range separation and using 
only the LDA for E£ c [n], the approach already yields 
good results for total energies of atomic and molecular 
systems [2l|, HBJ . The purpose of this paper is to ana- 
lyzed in further details the two separate needed approx- 
imations of the method: the wave function calculation 
of Eq. (JSJ) and the short-range exchange-correlation func- 
tional E% c [n] . Accurate calculations have been performed 
on a few small atomic systems to assess the approxima- 
tions. Concerning the wave function part, we show that 
the modification of the interaction enables to decrease 
the effort, or alternatively increase the accuracy, of a CI- 
type calculation. Concerning the functional part, we gain 
more insights into the short-range functional and explain 
the performance of the LDA by studying the behavior 
of E£ c [n] with respect to the interaction parameter /i. 
In particular, we show that, when the interaction v£ e is 
short ranged enough, E£ c [n] can be expressed as a sim- 
ple local functional of the on-top pair density and its first 
derivatives, explaining the accuracy of (semi)local func- 
tional approximations. This analysis also enables us to 
propose new approximations for Eg c [n] which correct the 
LDA in the domain of fi where it fails. 

The paper is organized as follows. In Sec.[nJ we give 
two possible choices for the long-range interaction v£ e (r) . 
In Sec. IIIII technical details concerning the calculations 
made on atomic systems are given. In Sec. lIVI the impact 
of the modification of the interaction on the performance 
of wave function calculation is investigated. In Sec. [V] 
we study the behavior of the corresponding short-range 
exchange and correlation functionals with respect to the 
interaction parameter /i. In Sec. lVll we test the obtained 
exact behaviors and compared them with the LDA. In 
Sec lVIII interpolations for the short-range exchange and 
correlation functionals with respect to \i are proposed 
and tested. Finally, Sec I Villi draws a conclusion of this 
work. 

Atomic units (a.u.) will be used throughout this work. 



J5" 



LDA r 



n(r)eZr il (n(r))dr, 



(9) 



II. LONG-RANGE /SHORT-RANGE 
SEPARATION 



where e^ u (n) is the complement short-range exchange- 
correlation energy per particle obtained by difference 
from the exchange-correlation energies per particle of the 
uniform electron gas with the standard Coulomb inter- 
action, e"' c llf (n), and with the long-range interaction vg e , 



€ 



/j,,unif 



(n) 



„unif 



(n) 



_/i,unif 



(n). 



(10) 



The method has been implemented at an experimen- 
tal level into the quantum chemistry package Molpro, 
allowing a combination of Cl-type wave function and 
DFT calculations |2l|. Recently, the method has also 



In a number of previous works [3 HHHI^ , a 
splitting of the Coulomb interaction based on the error 
function has been studied to describe the long-range part 
of the electron-electron interaction 



<e,erfO) 



erf(/ir) 



(11) 



This interaction, referred to as the erf interaction, has 
also been used by other authors in DFT for various pur- 
poses mmiMEHIillilHIi!. In this work > we 

introduce another interaction achieving a sharper sepa- 
ration of long-range and short-range interactions by sub- 
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tracting a Gaussian function from the erf interaction 



r V7T 



ee,erte;au 



where the coefficient and the exponent of the Gaussian 
have been chosen so that v£ e erfgau (r) and its derivative 
with respect to r vanish at r = 0. This modified in- 
teraction is referred to as the erfgau interaction. No- 
tice that this partition of Coulomb interaction has al- 
ready been proposed by Gill and Adamson 0. In 
another context, Prendergast et al. |3(| have also used 
a similar form of long-range interaction. For the sake 
of completeness, we finally note that another possible 
form of modified electron-electron interaction based on 
the Yukawa potential has also been investigated in the 

past 0111 m 13 mm. 

Both interactions l|ll|) and l|12ll enable to define a 
generalized adiabatic connection |39j| between the non- 
interacting KS system at fj, = where the interac- 
tion vanishes v^°(r) = 0, and the physical system at 
fi — > oo where the full Coulomb interaction is recovered 

«£r°°(r) = l/r. 

Notice that, although we have chosen the same no- 
tation for convenience, the parameter (i in Eq. is 
a priori independent of that of Eq. (|12[1 . Actually, in 
all the plots of the paper involving the erfgau interac- 
tion, we apply a scale factor to the interaction parameter 
of the erfgau interaction: /i — > c/i. The constant c is 
chosen so as to have the same leading term in the dis- 
tributional asymptotic expansion of the two interactions 
when fi — > oo (see Appendix Q which leads to the value 

1 /2 

c = (l + 6\/3) ~ 3.375. This also insures that the 
leading term in the asymptotic expansion of short-range 
Hartree, exchange and correlation energies for large p, is 
the same for the two interactions. In Fig. the erf and 
erfgau interactions are compared along with the Coulomb 
interaction. The scale factor on the parameter of the er- 
fgau interaction enables to define a common "cut-off ra- 
dius" giving the range of the interaction and defined by 
the inverse of the interaction parameter r c w l//i. For 
interelectronic distances larger than r c , the two modi- 
fied interactions reproduce the long-range Coulomb tail. 
Notice that short-range interactions are better removed 
with the erfgau interaction. 

We finally note that, for Gaussian basis set calcula- 
tions, the evaluation of the two-electron integrals cor- 
responding to the erf or erfgau interaction requires only 
simple modifications of standard algorithms for Coulomb 
integrals (see Appendix |A"|) . 



III. DETAILS ON ACCURATE CALCULATIONS 

We explain rapidly how the accurate data for atomic 
systems presented in this work have been obtained. 

For each system, accurate potentials w M (r) needed 
in Eq. JSJ and corresponding value of the universal 



v£ e (r) 
(a.u.) 




FIG. 1: Coulomb interaction 1/r (dotted curve), erf interac- 
tion (dashed curve) and erfgau interaction (solid curve) for 
fi = 1. In order to compare the two modified interactions, 
a scale factor is applied on the interaction parameter of the 
erfgau interaction: fi — > cfi with c = 3.375 (see text). It is 
possible to define a common cut-off radius r c ~ 1/fJ,. 



long-range functional F M are computed using the Licb's 
Legendre transform formulation of the universal func- 
tional E3 



F"[n] = maxjE^] - J n(r)^(r)rfr|, 



(13) 



where [v^\ is the ground-state energy of the Hamilto- 

nian T+V£ + ^2i v^^i)- If n i s chosen to be the physical 
density of the system, the maximum is reached for the 
desired and F 11 . In practice, an accurate density n 
is computed by multi-reference CI with single and dou- 
ble excitations (MRCISD) [4 lL |42( and the potential to 
optimize v^(r) is expanded as 



n 
i=l 



C 

5 

r 



(14) 



where Cj are the optimized coefficients, pi are some 
fixed integers (—1 or 2), 7^ are fixed exponents cho- 
sen so as to form an even-tempered basis set (typically, 
7i G [10~ 3 , 5.10 4 ]), and C is a constant enforces the cor- 
rect asymptotic behavior for r — > 00. For the Kohn-Sham 
case (n = 0), C = —Z + N — l whereas for finite 
C = — Z (N and Z are the electron number and nuclear 
charge, respectively). 

The maximization of Eq. (|13|l is carried out with the 
Simplex method 01 . For a given potential, E^lv^] is 
computed at MRCISD level using Molpro 44] with mod- 
ified two-electron integrals (see Appendix |A"|> . Beside the 
asymptotic behavior for r — * 00, v^(r) ~ C/r, the behav- 
ior of the potential at the nucleus r = 0, v^(r) ~ —Z/r, is 
also imposed during the optimization. Large one-electron 
even-tempered Gaussian basis sets are used for all sys- 
tems (see Refs. and l26l for more details). 

The standard universal functional F and the KS poten- 
tial vks (and thus the KS determinant $) are obtained 
as a special case for /1 = 0. The complement short-range 
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functional = F — F M can then be deduced and the 
short-range exchange and correlation energies, Eg and 
Eg, are obtained from Eqs. 10 and ©■ 

IV. PERFORMANCE OF THE 
APPROXIMATIONS FOR THE WAVE 
FUNCTION CALCULATION 

We now investigate the effect of the modification of the 
electron-electron interaction on the wave function part of 
the calculation. For this purpose, we evaluate the effi- 
ciency of approximate resolutions of Eq. © as follows. 

We first construct, for each /i, the Hamiltonian us- 
ing an accurate potential w M and compute accurately its 
ground-state energy, E" = at the MRCISD 

level. We then compute various approximate ground- 
state energies, Eg = ( 1 i<g\H fJ -\'i'g), by using approximate 
Cl-type wave functions tyg expanded into linear combi- 
nations of all the few Slater determinants generated from 
small orbital spaces S. The orbitals used are the natural 
orbitals of the Coulombic system calculated at the MR- 
CISD level. The accuracy of the approximation for tyg 
can be assessed by looking at the difference between Eg 
and E^ 

AE^. =Eg- E». (15) 

The differences AEg are plotted along the erf and erf- 
gau adiabatic connections in Fig. [21 for the He atom with 
the orbital spaces S = Is, S = ls2s and S = Is2s2p. 
One sees that, in the Coulombic limit /i — + oo, the re- 
duction of the orbital space leads to important errors in 
the energy. When [i is decreased, i.e. when the interac- 
tion is reduced, the errors due to limited orbital spaces 
get smaller and smaller. For instance, at /i = 1, using 
only the single-determinant wave function ^i s , leads to 
an error AE^ S of less than 0.005 Hartree. The erfgau 
interaction generally gives smaller errors than the erf in- 
teraction, except near fj, = where the non-monotonicity 
of the erfgau interaction leads to peculiar behaviors [24| . 
Anyway, we will not use the erfgau interaction in this 
region of very small fi. 

The case of the Be atom with the orbital spaces S = 
ls2s and S = Is2s2p is reported in Fig. |3 Because of 
the near-degeneracy of the 2s and 2p levels, the inclusion 
of 2p configurations in the wave function is important, 
quite independently of the electron-electron interaction. 
Indeed, the difference E^ s2s — E^ s2s2p remains large for al- 
most all /i's. On the contrary, the error of the calculation 
where the 2p orbitals are included, AE^ s2s2p , quickly falls 
off when fi is decreased. Again, for fi = 1 for instance, 
the error AE^ s2s2p given by the few-determinant Cl-type 
wave function ^i s2s2p is less than 0.005 Hartree. 

Therefore, the modification of the interaction enables 
to increase the accuracy of Cl-type wave function cal- 
culations, or equivalently for a fixed target accuracy, de- 
crease the effort of the calculation by reducing the orbital 
space. The crucial point for this effect to appear is the 



0.04 



(a.u.) 




U (a.u.) 

FIG. 2: Ground-state energy differences A_Eg = 
- where * M is an accurate wave 

function and are approximate wave functions generated 
from small orbital spaces S = Is, S = ls2s and S — Is2s2p 
along the erf (full curves) and erfgau (dashed curves) adia- 
batic connections for He. 




V (a.u.) 

FIG. 3: Ground-state energy differences A_Bg = 
l$l%\k»\^%) - where *" is an accurate wave 

function and ^% are approximate wave functions generated 
from small orbital spaces S = ls2s and S — Is2s2p along the 
erf (full curves) and erfgau (dashed curves) adiabatic connec- 
tions for Be. 

reduction of the electron-electron interaction compared 
to the Coulomb interaction and not really the long-range 
character of the modified interaction. 



V. BEHAVIOR OF SHORT-RANGE 
EXCHANGE AND CORRELATION 
FUNCTIONALS WITH RESPECT TO THE 
INTERACTION PARAMETER 

At /i = 0, the short-range exchange-correlation func- 
tional reduces to the standard exchange-correlation func- 
tional of the KS scheme, Eg=° = E xc , and in the limit 
/i — > oo, the short-range functional vanishes, Eg^*°° — 0. 
Near these two limits, the study of the behavior of Eg 
and Eg with respect to fi constitutes an analysis of the 
exchange and correlation functionals in term of the range 
of the interaction. Indeed, the behavior at fi — > tells 
us how the KS exchange-correlation functional responds 
when very long-range interactions are removed from it, 
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while the asymptotic expansion for fi — > oo gives the 
exchange-correlation functional associated to very short- 
range interactions. To release the text from mathemati- 
cal details, the full derivation of the expansions are given 
in the Appendices. 



A. Exchange functional for small /i 



In Appendix^] we show that the short-range exchange 
energy has the following formal expansion around /i = 



-. oo 



7r L — ' n! 

n=0 

2», 



x n 2 , x (r 1 ,r 2 )r^dr 1 dr2, (16) 



where E x is the usual KS exchange energy, n 2 , x (ri,r 2 ) 
is the exchange contribution to the pair density, and a n 
are coefficients depending of the interaction chosen and 
defined after Eq. I|B2|) . More specifically, for the erf in- 
teraction, this expansion writes 



n 2 , x {r 1 .r2)rl 2 dr 1 dr 2 - 



(17) 

where the term linear in /i comes from the normalization 
of the exchange hole [3jJ . For the erfgau interaction, the 
expansion is 



E^ f = E x 

x.crtgau x 



2^ 5 
450? 



n 2tX {r 1 ,Y 2 )r\ 2 dv 1 dv 2 



(18) 

The exchange energy with the erf interaction varies lin- 
early in /i near the Kohn-Sham end of the adiabatic con- 
nection (see Fig. On the contrary, the exchange en- 
ergy with the erfgau interaction varies very slowly like /i 5 
which implies that the corresponding curve is very flat 
near /j, = (see Fig. - The latter behavior is not unex- 
pected for a good long-range/short-range separation. In- 
deed, with the erfgau interaction, when /i increases near 
fi = only very long-range interaction effects are re- 
moved from the functional; they practically do not later 
the exchange energy of a finite system. With the erf 
interaction, the long-range/short-range separation is im- 
perfect and thus the short-range part of the interaction 
is also affected near /i = which is responsible for the 
linear behavior of the exchange energy in this case. This 
clearly shows that the erfgau interaction realizes a bet- 
ter separation of long-range and short-range interactions 
than the erf interaction. 



B. Correlation functional for small /i 

The general expansion of short-range correlation en- 
ergy around fi = is derived in Appendix [BJ It reads 



oo oo 




(-l)"(2n + l)a„ j„ +t+ i 

H 1 



7T ^ n\k\(2n + k + 1) 

n=l k=l v ' 



d k nt 2 l c (r 1 ,r 2 ) 
dfi k 



r^dridr 2 . (19) 



M=0 



where E c is the usual correlation energy of the KS scheme 
and n 2 c {yi,v 2 ) is the correlation pair density with in- 
teraction v£ e . Actually, several terms of this expansion 
vanish. For the erf interaction, the expansion writes 



Kerf = E c + 




360? 

dfi 3 



r? 2 dridr 2 + 



(20) 



Similarly, the expansion for the erfgau interaction is 
u 10 

E^ F — 

Ccrfgau c 5400 ^: 




d 5 < e (ri,r 2 ) . , 



r 12 dr\dv 2 + 



(21) 



The correlation energy varies much more slowly than the 
exchange energy near the KS end of the adiabatic con- 
nection. Therefore, the curve of the correlation energy 
with respect to fj, is very flat around /i = (see Figs. El 
and EJ). Again, this is due to the fact that removing very 
long-range interactions from the functional has no effect 
in a finite system. It will be seen in Sec IVHI that this 
makes the correlation energy difficult to interpolate near 
H = from a knowledge of the functional for large /i. 



C. Exchange functional for large /j 

In Appendix [CI we derive the general asymptotic ex- 
pansion of the short-range exchange energy for [i — > oo 



E£ = 2^£ 



A 2n 



(2n)!(2n + 2) M 



2n+2 



/n 2 2 : ) (r,r)dr + 0(-^), (22) 



where (r, r) arc the on-top exchange pair density 
and its spherical-averaged (with respect to r 12 ) deriva- 
tives. Simple explicit expressions can be given for the 
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first two terms of this expansion. Indeed, the on-top ex- 
change pair density writes 

n ailB (r,r)=-5^nJ(r), (23) 

where the summation is over the two spin states a = a, (3 
and n a (r) are the spin-densities; its second-derivative can 
be expressed by jig 

x (; Mr) _ 4Mr)+ M), 

(24) 

where r CT (r) are the KS spin kinetic energy densities ex- 
pressed in term of the KS spin-orbitals (r) by 

, N a 

^(r) = 2£l V< M r )! 2 - (25) 

i=l 

where N a is the number of electrons of spin a. The Lapla- 
cian in Eq. (|24|l can be eliminated by integration by parts 

J n a (r)V 2 n a (r)dr = - J \Vn a (r)\ 2 dr, (26) 

which leads for the leading terms of the expansion when 
fi — ► oo of the short-range exchange energy 

+ •••, (27) 

With A 0>e rf = y/n/2, A , C rfgau = (l + 6\/3)A),crf, A 2:Cr i = 

30F/4 and A 2 , C rfgau = (1 + 36v / 3)^2,Grf ■ Note that 
this last result has already been demonstrated for the 
case of the erf interaction |29| . In particular, for spin- 
unpolarized systems, Eq. (|27|l becomes 

V^4 2 f , , /|Vn(r)| 2 „ , A , 

(28) 

Eq. H27fl or 128|) shows that the first term in the expan- 
sion for large /i of the short-range exchange energy is an 
exact local functional of the density (or, in general, of 
the spin-densities). Therefore, the local density approxi- 
mation becomes exact in this limit. An alternative view 
of this result can be achieved in the framework of the 
wave- vector analysis of the usual KS exchange functional 



by considering the short wavelength limit [5|. The next 
term of expansion H28|l involves the gradient of the den- 
sity |Vn| and the kinetic energy density r and is therefore 
of the meta-GGA type. The following higher-order terms 
involve of course more and more ingredients constructed 
from higher-order derivatives of the KS orbitals. Note 
that instead of considering these spin-dependent quanti- 
ties one can also directly use the on-top exchange pair 
density and its derivatives. 

A similar asymptotic expansion can be derived for the 
short-range Hartree energy 




By carefully comparing Eq. (|27(l and Eq. (|29(l . one sees 
that, for a one-electron system (n a = n, np = 0, r Q = 
i~w = |Vn| 2 /(8n), Tp = 0), the Hartree and exchange 
energies cancel out order by order in the expansion with 
respect to fi. As the first term (in n~ 2 ) of the asymptotic 
expansion of the exchange energy is a local functional of 
the density, this cancellation is exactly maintained in the 
LDA, i.e. there is no self-interaction error in LDA for 
this first term. For the next term of the expansion (in 
/i~ 4 ), the removal of the self- interaction error requires 
the consideration of a meta-GGA functional depending 
explicitly on |Vn CT | and r a . 

In the KS scheme, it has been often argued that 
(semi)local approximate exchange functionals actually 
mimics near-degeneracy correlation effects (see, e.g., 
Ref. It?]) . It is clear from this above discussion that, 
for short-range exchange functionals, (semi)local approx- 
imations become exact at large fi and therefore do not 
fortuitously mimic near-degeneracy correlation anymore. 

D. Correlation functional for large /i 

In Appendix [CJ we derive the first two terms in the 
asymptotic expansion of the short-range correlation en- 
ergy for fi — > 00 

Eg = ^ 2 J n 2 ,c(r,r)dr+ ^ J n 2 (r,r)drH , 

(30) 

where «2,c(r, r) and n 2 (r, r) are the correlation and total 
pair density, the A coefficient is given after Eq. (|27|l and 
in addition Ai jCr f = 1 and ^l^rfgau = 28. 

This result show that for large /i the short-range cor- 
relation energy becomes an exact local functional of 
the on-top pair density n 2 (r,r) (and of the density via 
7i2, c (r, r) = n 2 (r, r) — n(r)n(r)/2 for closed-shell sys- 
tems) . This emphasizes the importance of this quantity 
for density functional approximations. It is in the same 
line of thought of a number of previous studies |43.l4^.l50| 
stressing the need of including explicitly the on-top pair 
density in approximate density functionals. 
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VI. PERFORMANCE OF THE LDA FOR 
SHORT-RANGE EXCHANGE AND 
CORRELATION ENERGIES 



The asymptotic expansions for /i — * oo enables to an- 
alyze the local density approximation to the short-range 
exchange-correlation functional. We restrict the discus- 
sion to spin-unpolarized systems. Applying the LDA to 
the expansion of the exchange energy (|22[) corresponds 
to transferring the on-top exchange pair density and its 
derivatives n^ 2 ™" 1 (r, r) from the uniform electron gas |5lT ] 



4!r ) ( r . r ) ~ n(r)n(r) (gf${n{r)) 



(31) 



where <^. 2 Q^(n(r)) are the on-top derivatives of the 
exchange-only pair-distribution function of the uniform 
electron gas £3 



9 /sin(fc F ri 2 ) - k F r 12 cos(k F r 12 ) 
9x{ri2 ' n) = 1 " 2 VWf, 



V F' 12 



(32) 

with k F = (37T 2 ?!) 1 / 3 . The expansion of the short- 
range exchange energy then writes 



EH 



LDA 



2p? 

3 2/ 37r ll/6 A2 

20/x 4 



nivfdv 



i(r) 8 / 3 dr 



(33) 



This last expansion is compared with the exact expan- 
sion (|28|l for the He atom with the erf and erfgau interac- 
tions in Figs. 01 and |5j respectively. In these expressions, 
an accurate density of the He atom is used. The LDA 
for the modified interactions (l7l l24j evaluated with the 
same density is also reported, as well as an accurate cal- 
culation of the exchange energy along the adiabatic con- 
nection [11113 (see Sec. HE). The LDA expansion Q 
is exact for the first term and very close to the exact ex- 
pansion l|28l) with the first two terms. Consequently, in 
the domain of validity of this expansion (for /i > 2), the 
LDA is nearly exact. The success of the LDA for large 
/i is therefore due to the exactness of the exchange on- 
top pair-density and the good transferability of its first 
derivatives from the uniform electron gas to the finite 
system. Actually, it can be remarked that the expansion 
of LDA up to the fi~ 4 term performs even slightly better 
than the exact expansion to the same order. Similarly, 
the LDA is already very accurate for /i > 1 while achiev- 
ing a comparable accuracy from the slowly-improving ex- 
pansion (|28|l of the exact exchange energy would require 
a rather long expansion. 

For the correlation energy, the LDA of expansion l|3U|) 
consists in transferring the on-top density from the uni- 



(a.u.) 




(a.u.) 



FIG. 4: Accurate short-range exchange energy of He (thick 
solid curve) along the erf adiabatic connection, local density 
approximation (thick long-dashed curve), exact asymptotic 
expansion for fj, — ► co (Eq. I28B with the first term (lower 
solid curve) and with the first two terms (upper solid curve), 
and asymptotic expansion in LDA (Eq. 1331 with the first two 
terms (short-dashed curve), the first term in LDA being exact. 



(a.u.) 




V (a.u.) 

FIG. 5: Accurate short-range exchange energy of He (thick 
solid curve) along the erfgau adiabatic connection, local den- 
sity approximation (thick long-dashed curve), exact asymp- 
totic expansion for /i — + oo fEq. 12811 with the first term (lower 
solid curve) and with the first two terms (upper solid curve), 
and asymptotic expansion in LDA (Eq. 1331 with the first two 
terms (short-dashed curve), the first term in LDA being exact. 



form electron gas to obtain 
sftA 



EH 



LDA 



^ 2 j n(r)n(r) lffo(r«(r)) - - j dr 
+ 2AiVtt ^ n ( r ) n ( r ) 5o ( rs ( r )) dr + ... ^ 



3a* 3 



(34) 



where r s — 3/(47m) 1 / 3 is the local Wigner-Seitz radius 
and go{f s ) is the on-top pair-distribution function of the 
electron gas for which Burke, Perdew and Ernzerhof have 
proposed an estimation |5J| 

g {r s )=D([ 1 + r s f/ 2 +p)e- A ^^, (35) 

with D = 32/(3tt), A = 3.2581, /? = 163.44 and 
7 = 4.7125. The expansion (|34|l is compared with the 
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FIG. 6: Accurate short-range correlation energy of He (thick 
solid curve) along the erf adiabatic connection, local density 
approximation (thick long-dashed curve), exact asymptotic 
expansion for fi — > oo (Eq. I3UI with the first term (lower 
solid curve) and with the first two terms (upper solid curve), 
and asymptotic expansion in LDA f Eg, 1341 with the first term 
(lower short-dashed curve) and with the first two terms (upper 
short-dashed curve). 
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fi (a.u.) 



FIG. 7: Accurate short-range correlation energy of He (thick 
solid curve) along the erfgau adiabatic connection, local den- 
sity approximation (thick long-dashed curve), exact asymp- 
totic expansion for /i — » oo (Ea. l3UH with the first term (lower 
solid curve) and with the first two terms (upper solid curve), 
and asymptotic expansion in LDA (Eg. 1341 with the first term 
(lower short-dashed curve) and with the first two terms (up- 
per short-dashed curve). 



exact expansion (|30|) using an accurate calculation of 
n.2 (r, r) for the He atom with the erf and erfgau inter- 
actions in Figs. EJand [7] respectively. An accurate calcu- 
lation of the correlation energy computed along the adi- 
abatic connection ptil l53| and the LDA for modified in- 
teractions 0, 13 are also reported. As for the exchange 
energy, the first terms of the LDA expansion (|34|) nearly 
coincide with the exact expansion. This expansion with 
the first two terms gives a very accurate approximation 
to the exact correlation energy in the region of accuracy 
of the full LDA curve (from \i ps 2 to — > oo). From 
these results, it is clear that the total on-top pair density 
«2 (r, r) have good transferability from the uniform elec- 
tron gas to the He atom. Actually, this good (while not 
exact) transferability seems quite general and has already 
been pointed out for several atomic and molecular sys- 
tems [5(J, Hj, |55| • This gives an explanation for the suc- 
cess of the LDA in treating short-range electron-electron 
interactions (large fi) . On the contrary, toward the Kohn- 
Sham end (/i = 0) of the adiabatic connection, the LDA 
transfers spurious long-range correlations from the uni- 
form electron gas and therefore poorly extrapolates the 
exact correlation energy of the finite system which does 
not contain these long-range correlation effects. 



beyond the LDA such as gradient-corrected functionals. 
A simple idea for improving the short-range exchange- 
correlation energy functional along the adiabatic connec- 
tion is therefore to interpolate between an available den- 
sity functional approximation (DFA) for /i = and the 
/i-dependent LDA for /i — > oo. In the spirit of the usual 
DFT approximations, this interpolation will be done lo- 
cally, i.e. for the short-range exchange-correlation energy 

n] via 



density £% c (r) related to the global functional 



drn(r) e£ c (r). 



(36) 



A. Rational interpolations 

A simple possibility is to interpolate e£ c (r) along the 
adiabatic connection using an estimate at \x = and 
the expansions for large (i (and eventually for small fi) 
presented in Sec. For example, consider the rational 
approximant for the short-range exchange energy density 
of the erf interaction 



-DFA 



VII. INTERPOLATIONS FOR THE 
SHORT-RANGE EXCHANGE AND 
CORRELATION FUNCTIONALS 

In the previous section, we have shown that the LDA 
treats successfully short-range interactions correspond- 
ing to large interaction parameters fi but is inaccurate 
toward the KS end of the adiabatic connection, i.e. for 
small fi. But for fi — a lot of better estimates of the 
exchange and correlation energies are available with den- 
sity functional approximations of the KS scheme which go 



1 + hfl + b 2 ^ 2 



(37) 



where h = -l/(e£ FA 0F) and b 2 = -Ae^/^n) are 
chosen to satisfy the expansion for small [i (Eq. 117(1 and 
the expansion for large /x (Eq. I28[) to leading order. In 
Eq. (33, e° FA is the exchange energy density at fj, = 
estimated by one of the usual density functional approx- 
imations of the KS scheme. 

For the erfgau interaction, the short-range exchange 
energy density can also be interpolated between e^ FA at 
H = (no linear term in //, Eq. I18|) and the expansion 



9 



for large n (Eq. I28|) according to 



-DFA 



""a;, erfgau 



1 



(38) 



with c 2 - -4e° FA /((l + 6V5)tto). 

In the same spirit, the short-range correlation density 
along the adiabatic connection can be approximated by 



-DFA 



1 + dill + d,2\i 



(39) 



with dx^xi = -2rf2,crfV / ^ rl 9o(^)/(3e c ), d 2 ,orf = 
2e c /(irn(g (r s ) — 1/2)) for the erf interaction, and 
dixrfgau = -56dl er{gan ^/Trng (r s )/(3e c ), <i 2 , erfgau = 
^2,crf/(l + 6v / 3) for the erfgau interaction. Again, eJ? FA is 
the correlation energy density at /i = given by an avail- 
able density functional approximation of the KS scheme. 
The rational approximant (|39|l therefore interpolates be- 
tween e^ FA at [i = and the expansion of the LDA for 
large n with the first two terms (Eq. I34fl . 



B. Weighted interpolations 

The interpolation formulas l|37l) and l)38|) only make 
use of the first term of the asymptotic expansion of the 
LDA exchange energy for /j, — * oo. However, it has been 
realized in Sec. IVII that the LDA works well on a larger 
domain of [i than the first terms of its asymptotic expan- 
sion. In order to take better advantage of the LDA, one 
can modify it only in the region of small /i by using infor- 
mation from better estimates of the KS exchange energy 
at /j, = 0. We therefore interpolate locally the short-range 
exchange density along the adiabatic connection by 



DFA 



£ 



p— 0,unif 



Mm) 



(40) 



where w([i) is a weight function acting at small fj, only. 
More precisely, w(/z) must be a positive function satisfy- 
ing w(fi = 0) = 1 so that e^ =0 = e^ FA , having significant 
values in the region where the LDA fails and with a fast 
decay for fj, — ► oo so as to recover the correct behavior 
of e^' unlf for large fi. For the erf interaction, we found 
that the local value of fi delimiting the domain where 
the LDA is inaccurate is well estimated by l/r s (r) where 
r s (r) is the local Wigner-Seitz radius. We therefore take 
w m f (/i) = erfc(/zr s ) where the complementary error func- 
tion erfc ensures that the weight function has significant 
values only for fi < l/r s . Eq. i|4U|) thus interpolates be- 
tween eJ? FA at m = and the LDA at large [i. 

For the erfgau interaction, we can used the same weight 
function except that now the local value of [i delimiting 
the domain where the LDA must be corrected is esti- 
mated by /i — c/r s where c = (l + 6\/3) ~ 3.375 is 
the scale factor between the erf and erfgau interactions 
discussed in Sec. [HI We therefore take as weight function 
w er fgau(^) = erfc(/ir s /c). 



Naturally, this interpolation can also be applied to the 
short-range correlation energy 



DFA 



-/j,— 0,unif 



(41) 



with the same weight function w(/i). 



C. Other interpolations 

Finally, we mention the approximation for the short- 
range exchange e nerg y proposed by Iikura, Tsuneda, 
Yanai and Hirao based on a modification of the 
short-range LDA exchange functional [TtI l24l | . The spin- 
unpolarized version of their approximation for the erf in- 
teraction is 



.DFA 



-(2A - 4A 3 )e~ 1/{4A2) -3A + 4A 3 



,(42) 



where A 



/i/(2fc), k = y/el; DA /e^ FA k F and k F 

-DFA 



(S^n) 1 / 3 . The approximation reduces to £^ rA at fi = 
and has an asymptotic expansion for /i — > oo incorpo- 
rating the correct leading term (cf. Eq. I28JI . There- 
fore, Eq. 1)42(1 provides an interpolation between a density 
functional approximation at fi = 0, £^ FA , and the correct 
limit as [i — ► oo. The same approximation can also be de- 
rived for the erfgau interaction using the LDA exchange 
functional associated to this interaction 241; it reads 



.DFA 



""x, erfgau 



+(2A - AA 3 ^ 1 '^ - 3A + AA 3 

+ {2B - 16B 3 )e- 1/(4s2) -6B + 16B 3 ^j 

(43) 

where B = A/^/3. 

N otice t hat in the interpolations of Sec. IVII Al IVII Bl 

and IVII O the exchange and correlation energy densities 
at = 0, eJ? FA and e° FA , can be estimated by any of the 
available exchange-correlation functionals of the Kohn- 
Sham scheme. In this sense, formulas (13711 to (|43|) provide 
extensions of these exchange-correlation functionals over 
the erf and erfgau adiabatic connections. 



D. Results 

We now test the interpolations formulas l|37|l to l|43|) 
on a few atomic systems. In all the results presented 
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FIG. 8: Accurate short-range exchange energy of He (solid 
curve) along the erf adiabatic connection, local density ap- 
proximation (long-dashed curve) and three interpolations be- 
tween the PBE value at fj, = and the LDA at fi — » oo 
using a rational approximant (Eg. 1371 short-dashed curve), a 
weighted approximant (Eq. l4UI dotted curve) and the approx- 
imation of Iikura el al. (Eq. 1421 dashed-dotted curve) . The 
last two curves are nearly superimposed. 
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FIG. 9: Accurate short-range exchange energy of He (solid 
curve) along the erfgau adiabatic connection, local density 
approximation (long-dashed curve) and three interpolations 
between the PBE value at fj, = and the LDA at y, — * oo 
using a rational approximant (Eg. 1381 short-dashed curve), a 
weighted approximant (Eo. l4L)l dotted curve) and the approx- 
imation of Iikura el al. (Eq. 1431 dashed-dotted curve) . The 
last two curves are nearly superimposed. 



here, we use the PBE functional j)6] as the density func- 
tional approximation at fi — both for exchange and 
correlation: £ ° FA = £ p be and e» FA = e? BE . 

Figs. |H1 and represent the short-range exchange en- 
ergy of the He atom along the erf and erfgau adiabatic 
connections, respectively. An accurate calculation is 
compared to the LDA and to the three interpolations 
of Sec. IV11AI IVlTBl and IVrTTl One sees that the ra- 
tional approximants give an overall reasonable estimate 
of the exchange energy along both adiabatic connections 
but are actually less accurate that the LDA in the re- 
gion of intermediate [i. As already noticed, this reflects 
the fact that the LDA works well on a larger range of 
/i that the first terms of its expansion for fi — > oo. The 
weighted interpolations make better use of the LDA and 
constitutes an improvement over it for all /x's. Finally, 
the approximation of Iikura et al. is nearly identical to 
the weighted interpolation for both the erf and erfgau 
interactions. 

Similar curves for the short-range correlation energy 
of the He atom with the erf and erfgau interactions are 
shown in Figs. El and ^2 As for the exchange energy, 
the rational approximants constitute an overall correc- 
tion to the LDA but not in a systematic way since the 
LDA still performs better for intermediate fi. On the 
contrary, the weighted interpolations always improve the 
LDA. As one can expect, considering exchange and cor- 
relation together further improves the results at small fi. 

We now discuss the Be atom. The short-range ex- 
change energy along the adiabatic connections looks very 
similar to that of the He atom and will not be shown. 
The short-range correlation energy for which the LDA 
has more difficulties than for the He case is represented 
in Figs. El an d El f° r the erf and erfgau interactions. 
One can see that both the rational approximant and the 
weighted interpolation improve the LDA along the whole 
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FIG. 10: Accurate short-range correlation energy of He (solid 
curve) along the erf adiabatic connection, local density ap- 
proximation (long-dashed curve) and two interpolations be- 
tween the PBE value at n — and the LDA at n — > oo using 
a rational approximant (Eq. 1391 short-dashed curve) and a 
weighted approximant fEo. 1411 dotted curve). 



adiabatic connection. 

The He and Be atoms are simple cases where a 
gradient-corrected functional like PBE give a very accu- 
rate correlation energy. On the contrary, the Ne 6+ atom 
constitutes a much more difficult system for (semi) local 
functionals because of the presence of strong near- 
degeneracy correlation effects due to the proximity of the 
2s and 2p levels. The short-range correlation energy of 
this system is reported in Figs. El an( i El f° r the erf 
and erfgau interactions. The PBE functional of the KS 
scheme (/i = 0) strongly underestimates the correlation 
energy. One sees that with the rational approximants or 
the weighted interpolations, the error is rapidly decreased 
when fi is increased, i.e. when long-range interactions are 
removed from the functional. 
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FIG. 11: Accurate short-range correlation energy of He (solid 
curve) along the erfgau adiabatic connection, local density ap- 
proximation (long-dashed curve) and two interpolations be- 
tween the PBE value at /i = and the LDA at fj, —* oo using 
a rational approximant (Eq. 1391 short-dashed curve) and a 
weighted approximant (Eq. ^3 dotted curve) . 
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FIG. 14: Accurate short-range correlation energy of Ne 6+ 
(solid curve) along the erf adiabatic connection, local density 
approximation (long-dashed curve) and two interpolations be- 
tween the PBE value at = and the LDA at n — > oo using 
a rational approximant (Eq. 1391 short-dashed curve) and a 
weighted approximant fEq. 1411 dotted curve). 
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FIG. 



12: Accurate short-range correlation energy of Be (solid 
curve) along the erf adiabatic connection, local density ap- 
proximation (long-dashed curve) and two interpolations be- 
tween the PBE value at ii — and the LDA at /i — > oo using 
a rational approximant (Eq. 1391 short-dashed curve) and a 
weighted approximant fEq. 1411 dotted curve). 
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FIG. 15: Accurate short-range correlation energy of Ne 6+ 
(solid curve) along the erfgau adiabatic connection, local den- 
sity approximation (long-dashed curve) and two interpola- 
tions between the PBE value at fi = and the LDA at fi — > oo 
using a rational approximant fEa. 1391 short-dashed curve) and 
a weighted approximant (Eq. 1411 dotted curve) . 
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FIG. 13: Accurate short-range correlation energy of Be (solid 
curve) along the erfgau adiabatic connection, local density ap- 
proximation (long-dashed curve) and two interpolations be- 
tween the PBE value at /i = and the LDA at fj, —* oo using 
a rational approximant (Eq. 1391 short-dashed curve) and a 
weighted approximant (Eq. 1411 dotted curve) . 



VIII. CONCLUSION 

In a quantum electronic system, the long-range/short- 
range separation of the Coulomb interaction enables to 
rigorously decompose the total energy into long-range 
and short-range components which can be calculated by 
different methods. In particular, a density functional ap- 
proximation can be used for the short-range part of the 
energy, while the long-range contribution can be treated 
by traditional wave function methods. In this work, we 
have considered two possible long-range/short-range sep- 
aration of the Coulomb interaction: the erf and erfgau 
modified interactions. The erfgau interaction achieves 
a better separation than the erf interaction. We have 
shown that the use of these modified interactions facil- 
itates wave function calculations. We have also studied 
the short-range part of the exchange-correlation func- 
tional with respect to the range of the associated inter- 
action, and shown that, in the limit of a very short-range 
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interaction, the exchange contribution to this functional 
can be expressed as a local functional of the density and is 
therefore exact in the LDA. In the same limit, the correla- 
tion contribution can be expressed as a local functional of 
the on-top pair density which is generally accurate in the 
LDA. However, when the interaction becomes more and 
more long-ranged, the LDA is inaccurate. It is neverthe- 
less possible to improve the LDA description of the short- 
range exchange-correlation functional by making use of 
the available gradient-corrected functional of the Kohn- 
Sham scheme in the limit of the full Coulomb interac- 
tion. We have indeed proposed two kind of interpola- 
tions for the short-range exchange-correlation functional 
along the erf or erfgau adiabatic connection which im- 



prove the LDA. In order to extend the interaction range 
well treated by functional approximations in a more sys- 
tematic way, we are currently investigating gradient cor- 
rections for the short-range exchange-correlation func- 
tional. 
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I 

APPENDIX A: MODIFIED TWO-ELECTRON INTEGRALS 



In this Appendix, we give details of the evaluation of the two-electron integrals over Gaussian basis sets for the 
modified erf and erfgau interactions. The modified integrals have been implemented into the Seward program [57j 
available in the Molpro package | l 1| and using the Rys quadrature scheme [58(. For other implementations of two- 
electron integrals with modified interactions, see Refs. l28l l36l \E$L 

The general form of four-center two-electron integrals is 

[ab|u ee (ri 2 )|cd] = J J dridr 2 (/>(r 1 , £ Q , a, A)0(ri, Cb, b, B)u e e(?"i 2 )0(r 2 , Co c, C)0(r 2 , Q, d, D), (Al) 

where v ee {ri2) is the electron-electron interaction and <j) is is an unnormalized primitive Cartesian Gaussian basis 
function 

0(r,C,n,R) = {x-R x ) n *{y-Ry) n y{z-R z ) n *e-^ r -^ 2 . (A2) 
Actually, it is only required to explicitly compute the integrals of the type 

[eOK(r 12 )|fO] = k ab k G d ([ dv 1 dv 2 {x 1 - A x ) e *{ Vl - A y ) e * (zi - A z f* e -^~ F ^ 'v ee {r 12 ) 

J J (A3) 

(x 2 - C x y*(y 2 - C y yy(z 2 - e-^-Q)", 

where C = Ca + Cb, V = Cc + Crf, P = (CaA + C & B)/(C a + Cb), Q = (CcC+CdD)/(Cc+Gi), Mb = exp{-[CaC&/(Ca+C&)](A- 
B) 2 } and kqd — expj— \Cr.(d/(Cr. + G)](C — B) 2 |. The general integrals [ab|w ee (ri 2 )|cd] are then obtained by applying 
the so-called transfer relation (53 to the integrals [e0|u ee (ri2)|f0] with e = max(a, b)...a+b and f = max(c, d)...c + d. 
For the Coulomb interaction v ee (ri 2 ) = V r i2> these integrals can be written as 

[eO|l/r 12 |fO] = f dtP n (t)e- t2 ^ p -^ 2 , (A4) 
Jo 

with p — Cv/(C + v) an< i Pn(t) are the Rys polynomials of order n = e x + e y + e z + f x + f y + f z in t 2 

1/2 / \ 'V^ / \ 3/2 

P n {t)=2[^) n AB Kc D {-A f~J I x (e x ,f x ,t)l' v (e y ,f y ,t)l' z (e z J z ,t) (A5) 
where I' x (A = x, y, z) are 2D integrals obeying the recurrence relation 

I' x (e x + l,f x ) = f(P x -A x ) + ^(Q x -P x )\l x (e x J x ) 

+ | (l - J-) ^ - X > M + ^^ (eA ' h ~ X) ' (A6) 
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with the starting value I{(0, 0) = 1. Equation (|A4|) can be evaluated exactly by the Rys quadrature 

1/2 / \ 3 / 2 / \ 3 / 2 nR y 

[eO|l/r 12 |fD] = 2 kabkcd (~) (j) I x (e x J x ,t a )I y (e y J y ,t a )I z (e z J z ,t a )w a , (A7) 

where t a and w a = exp{— £ 2 p(P — Q) 2 } are the roots and the weights of Rys polynomials and nn ys > n/2. 

The modified integrals for the erf interaction [e0| erff/iri^/r^lfO! can simply be computed with the same scheme 
by applying the following simple modification everywhere [13, |2l||: l/p 1/p + 1/a* 2 - Thus, with respect to the 
Coulomb case, in Eq. I|A6|I the recurrence coefficients are modified, and in Eq. (|A7|I the prefactor, the weights w a and 
the roots t a , depending on p, are modified. 

To calculate the integrals for the erfgau interaction, additional integrals of type [eO|Cexp(— ar 2 2 )|fO] are needed. 
These integrals can be computed like an one-point Rys quadrature with root t a — \/a 2 / (p + a 2 ), weight w a — 
exp{— t 2 p(P — Q) 2 } and a modified prefactor depending on t a 

[eO|Cexp(-ar 2 2 )|fO] = C{1 - t 2 a f' 2 k AB k C d 

xl' x (ex, f x ,t a )I'(e y , f y , t a )I' z (e z ,f z ,t a )w a . (A8) 




APPENDIX B: ADIABATIC CONNECTION 
NEAR THE KS SYSTEM 

In this appendix, we study the erf and erfgau adiabatic 
connections near the KS system, i.e. for small interaction 
parameter p. The fictitious system along these connec- 
tions is described by the Hamiltonian [cf. Eq. (JSJ)] 

A " = E + X>-(^-) + £>"(*), (Bl) 

i i<j i 

where vg e is the long-range erf or erfgau interaction and 
is the external local potential associated to the short- 
range Hartree, exchange and correlation energy func- 
tionals v»(r) = v ne (r) + 5U»[n]/8n{v) + E£[n]/6n(r) + 
Eg[n]/6n(r). The behavior for // — > of all quantities 
associated to this fictitious system can be derived from 
the Maclaurin series of the long-range interaction vg e 



_ ^ ( 1 )" a " „2n,,2rt+l 
— > '12 P 



(B2) 



where a„ !0I f = l/(2n + 1) for the erf interaction and 
a n , erfgau = l/(2n + 1) — 1/3" ( = for n < 1 by con- 
struction) for the erfgau interaction. 



where U is the standard KS Hartree energy. To obtain 
Eq. (|B3() . the integral and summation signs have been 
interchanged. This is reasonable for a finite system where 
r*i2 is always bounded since in this case the series <|B2() 
is uniformly convergent and can thus be integrated term 
by term. 



2. Short-range exchange energy for p — ► 

The expansion for p — > of the short-range exchange 
energy of Eq. (7J is 



1 ( l) " a n 2n+l 



7r * — ' n! 

n=0 



x / / n 2 , x (Y 1 ,Y 2 )r{ 2 1 dv 1 dv 2l 



(B4) 



where E x is the standard KS exchange energy and 
^2,x( r ii r 2) is the exchange contribution to the pair den- 
sity. 



1. Short-range Hartree energy for yu — >■ 



3. Short-range correlation energy for p — > 



Eq. (|B2|1 leads immediately to the expansion of the 
short-range Hartree energy 

U' 1 = \JJ n{Y 1 )n{r 2 )v^ e {r 12 )dY 1 dv2 
_ tj 1 2n+l 



fix * — ' n! 

n=0 



x / / n(ri)n(r 2 )r 2 2 l dridr 2 



(B3) 



The short-range correlation energy of Eq. (JSJ can be 
written with the correlation contribution to the pair den- 
sity ri2 C (ri,r2) with interaction vg e as 

K = E --\\\ <c('i. r 2 )< e (r 1 2)dr 1 dr 2 , (B5) 

where E c is the usual correlation energy of the KS 
scheme. The expansion of Eff for small p can be ob- 
tained by first studying the derivative of Eg with respect 
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to fi which by the Hellmann-Feynman theorem writes 



dE? 1 



dfi 2 



< c (r 1; r 2 )^||pl dridr2i (B6) 



and assuming for n.2 C (ri,r2) the following expansion 
around /! = 



the correlation contribution. The two superscripts in 
^(2n+i,fc) re g ec ^ £ ac ^. ^jj a ^. terms in the expan- 



sion of the correlation energy (Eq. IB9(1 come from two 
sources: the interaction v£ e and the pair density n 2c . 
The ground-state wave function of this Hamiltonian 
can also be expanded with respect to fi 



~ 1 ( d k n^ c (r 1: r 2 ) \ k 

<c(n,r 2 ) = Efc! ^ ^ (B7) 

fe=i \ p / p=0 

since for the KS system n27°( r i' r 2) — 0. Therefore, 
inserting equation l|B7l) and the derivative of Eq. I|B2|) 



into l|B6(l , and assuming the commutativity of summation 
and integration, leads to 



dEg 



^ oo oo 



(-l )"(2n + l)o n 2n+fe 
/7r 4 — ' ' — ' n\k\ 

n=0 k=l 




d*n£ e (n,r 2 ) 
d[x k 



r^dridr 2 , 



/i=0 



(B8) 



where the term n = can be dropped since 7i2 C (ri,r2) 
integrates to zero. After integration we obtain the ex- 
pansion of the short-range correlation energy 

™ _ F 1 ^ y> (-l) n (2n + IK 2 »+fc+i 
" ' c 7i!fc!(2n + fc+ 1) M 




9 fc n^ e (ri,r 2 ) 



r^dndr 2 . (B9) 



M=0 



Actually several terms of this expansion vanish as it will 
be seen below. 



4. Ground-state wave function for /i — * 

Using the expansion of the interaction t^ e ('Ea. lB2l) and 
taking the functional derivatives of the expansions of the 
short-range Hartree, exchange and correlation energies 
(Eqs. IB3I IB4I and IB9|) to obtain the expansion of the 
potential v M , we arrive to the following formal expansion 
for the Hamiltonian of the fictitious system 

oo oo 
H" = J ffA-5 + E^e n+1) M 2 " +1 + E^ n+1) M 2 " +1 



n=0 
oo oo 



71=0 



+ l+k 



n=l k=l 



(BIO) 



(2n+l) 



In this equation, Hks is the KS Hamiltonian, v, 
refers to the coefficients of expansion (IB21. v {2 " +1) refers 
to the Hartree and exchange contributions to the coef- 
ficients of the expansion of and Vc +1 refers to 



fe=i 



(fe),,fc 



(Bll) 



where $ is the KS determinant. We now show that sev- 
eral terms in the expansion of actually vanish. 

For the erf interaction, vie* and are constants (cf. 

Eas. lBH lB3l andlBl | . and consequently v c l ' k) = for 
all k > 1, thus = 0. As there is no term in /i 2 



in i|B10(l . we also have \E'( 2 ) = 0. It implies in turn that 
in Eq. (|B7|I the terms k < 2 vanish, and thus by Eq. I|B9|) 



Vc +1 — for k < 2 and all n > 1. In particular, 
Vc 3 ' 1 ^ = Vc 3 '' 2 ^ = which leads to the expansion of the 
Hamiltonian 

K { = Hks + («2> + «& } ) + («<? + 

+w( 3 ' 3 V 6 + • • • . (B12) 

We therefore have the following expansion for the wave 
function 

*£ rf = $ + M 3 *( 3 ) + M 5 1r( 5 ) + M 6 *( 6 ) + • • • . (B13) 

Therefore, the terms k = 1,2,4 in Eqs. (|B"7jl and jB^ 

vanish for for the erf interaction. 

For the erfgau interaction, v^ee — vie — 0, and con- 
sequently = = and vi 3 ' k ^ = for all k > 1, 
thus *W = *( 2 ) = *( 3 ) = *( 4 ^ = 0. It implies in turn 
that in Eq. (|B7(1 the terms k < 4 vanish, and thus by 
Eq. (EH ^ 2 ™ +1 ' fe) = for fc < 4 and all n > 1. In par- 
ticular, Dc 5 ' 1 "* = Wc 5 ' 2 '' = Vc 5 ' 3 '' = fi 5 ' 4 " 1 = Vc = Vq 6 '' 2 ^ = 

(6,3) (6,4) (7,1) (7,2) (7,3) (8,1) 

Vc = V c = Vc — Vc — Vc = Vc = 

vi 8 ' 2 * 1 — Wc 9,1 ' ) = which leads to the expansion of the 
Hamiltonian 



^erfgau — H 



+ (^ ! +^)A 9 +W + '-. (B14) 

We therefore have the following expansion for the wave 
function 



crtgau 



V°* (10) + ' ' • ■ (B15) 



Therefore, the terms k = 1,2,3,4,6,8 in Eqs. jHJJl 
and l|B9|l vanish for the erfgau interaction. 



15 



APPENDIX C: ADIABATIC CONNECTION 
NEAR THE PHYSICAL SYSTEM 

In this appendix, we study the erf and erfgau adia- 
batic connections near the physical system, i.e. for large 
interaction parameter /x. 



1. Short-range interaction for fi — > oo 

We start by deriving a distributional asymptotic ex- 
pansion for large \i of the short-range electron-electron 
interaction vg e (r) = 1/r — v£ e (r) for the erf and erfgau 
interactions, 



c (r) 

ee,crfV / 



erfc(/ir) 



^ee, erfgau v ) 



eric(pr) 2^ i 2 r 



(CI) 



(C2) 



Let / : M + -> R be a test function (i.e., of bounded sup- 
port and infinitely differentiable) and consider the follow- 
ing integral 



1 = J /(r)<(r)dr. (C3) 
/ can be expanded into its Maclaurin series 

(C4) 



™ f(«)(0) 

' T7.I 



n=0 



where the Lagrange remainder R m (r) is 

(m + iy. 



(C5) 



with < 6 < 1. Inserting the first term of the right- 
hand-side of expansion 1|C4|> into l|C3(l gives the first con- 
tribution to integral / 



71=0 



(C6) 



where the last integral can be easily evaluated 



r n v% e (r)dr = 4tt / r n+2 v£ e (r)dr 



4^ 



(n + 2) /x"+ 2 



(C7) 

with ^n.orf = IX 1 ^) for the erf interaction and 
A^crfgau = r(2±5) + 3^(n + 2)r(2±2) for the erfgau 



interaction. Thus, Eq. (|C6H becomes 



n=0 



n\(n + 2) \i n 



+2 1 



(C8) 



The second contribution to / coming from the remainder 
writes 



h = 



(m + 1)! 



(C9) 



We shall assume in addition that, for any m, the (m + 1)- 
th derivative of / is bounded, i.e. |/( m+1 ) (r)| < M m , 
then 



\h I < 



r m+1 v£ e (r)dr 



(m + l)\ 

rn-\-l 

(m + l)!(m + 3) ^ m + 3 ' 



(CIO) 



meaning that I2 = C(l//i m+3 ). Finally, using the defini- 
tion of the n-th derivative of the three-dimensional Dirac 
delta distribution S^ n '(r), 



/(r),5(")(r)rfr=(-ir/ ( " ) (0), (Cll) 



we obtain the following distributional asymptotic expan- 
sion of v£ e (r) when fi — > 00 

<(r) = 4v^V ,/ } , ,(n) (r) (_J_). 



^ n!(n + 2)/i" 



(C12) 

Note that if we apply a scale factor to the interaction 
parameter of the erfgau interaction fi — ► cfx such that c = 

V^n.erfgauMn.erf) = (1 + 6^3) 1/2 w 3.375, then the 
asymptotic expansions of the erf and erfgau interactions 
are the same to leading order. This provides a criteria 
for comparison of the two interactions. 

Following the same procedure on the derivative of the 
modified interaction with respect to [i which writes 



d^erffr) 



2 22 



(C13) 



for the erf interaction, and 



^ee.erfgauM 



-fj, r 



dfx 



-.e 3' 



7T 

2„2 



4/rr 



e 3' 



(C14) 



for the erfgau interaction, leads to a similar asymptotic 
expansion 



(C15) 

which is just the derivative of expansion ilC12() . In the 
following, we will apply Eqs. {CTI and (|Cl5l) assuming 
that the corresponding test function /(r) related to the 
pair density satisfies all the required assumptions which 
are reasonable for finite systems. 
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2. Short-range Hartree energy for fi — » oo 

Introducing the trivial variable transformation r 2 — > 
r 12 and the spherical average of the density n(r 12 ) 

Hru) = ^ J n(r 12 )cM ri2 , (C16) 

the short-range Hartree energy writes 

m = \jj "( r i)"( r i2X( r i2)*idr 12 . (C17) 

Using then the distributional asymptotic expansion when 
fi — ► 00 of the short-range interaction (Eq. IC121 and 
noting that n(ri 2 ) can be expanded as an even series of 
r\2 around 7*12 = 0, we obtain the asymptotic series of 
the short-range Hartree energy 



£/" = 20F 



A 2n 



^ (2n)!(2n + 2)/i 2 «+ 2 



x / n(r)n^(r)dr + 0(^— (C18) 



1 



2m+3 - 



where the notation n^ 2n \r) — {d 2n h{r 12) / dr 2 2 l )r 12 =o 
has been used for the density and its spherical-averaged 
derivatives. 



3. Short-range exchange energy for /1 — > 00 

Similarly, the short-range exchange energy can be writ- 
ten as 

E *=\JJ "2,*(ri, ri 2 )i£ e (ri 2 )dridri 2 , (C19) 

where n 2;2; (ri, ri 2 ) is the spherical-average exchange pair 
density which can be expanded like the density as an even 
series of ri 2 around r\i = ("no cusp for exchange"). We 
therefore obtain the asymptotic series of the short-range 
exchange energy 



K = 2V^ 



A 2n 



x fn 2 2 :\r,r)dr + 0(^^), (C20) 
J A* 

with the on-top exchange pair density and 

its spherical-averaged derivatives n 2 2 ^ (r, r) = 

(0 a "na 1 *(r,n 2 )/0»ia) ri3 =o. 

4. Short-range correlation energy for n — > 00 

The asymptotic expansion of the short-range correla- 
tion energy when /i — > 00 can formally be found by con- 
sidering its derivative (cf. Eq. IB 61) 



dE{f 1 



dfi 



"2, C ( r l: r 12) 



dr lC ir 12 , (C21) 



where h 2 c (ri, r i2 ) is the spherical average of n 2 c (ri, r 2 ), 
inserting the asymptotic expansion of c (ri., ri.2) an d 
the distributional asymptotic expansion of dv^ e /dfi 
fEa. lC15j) and re- integrating with respect to fi. Actually 
considering only the first term of the asymptotic expan- 
sion of the correlation pair density i.e. ri 2c (ri,ri 2 ) = 
^2,c( r ii n.2) + " • • i is sufficient to find the first term of 
the asymptotic expansion of the correlation energy. Sim- 
ilarly to the exchange energy, the leading term of the 
short-range correlation energy for large \x is consequently 
given by the on-top correlation pair density n 2jC (r, r) 



= ^ I „ L ,,. (1 ..rWr- 



(C22) 



Notice that, for the erf interaction, this result has already 
been derived |2r3 |. 

Using Eqs. ifCTS)) . ltC20|) and l(C22|l . the external local 
potential u M (r) = v ne (r) + SU^ /8n(r) + 5E£ c /5n{r) has 
the following behavior for large p 



v tl (r)=v ne (r) + 



J ra 2 (r',r')dr / + 



[i 2 Sn(r) 

(C23) 

where n 2 (r, r) is the total on-top pair density. As v£ e has 
a similar expansion for /1 — > 00 beginning with pT 2 (see 
Eq.EEl) 



r 



(C24) 



the behavior for large \i of the modified Hamiltonian 
H^ = f+ v&fa) + v^Ti) is of the form 



H" =H+—H {2) 
H 2 



(C25) 



The corresponding wave function and consequently 
the correlation pair density n 2c (ri,ri 2 ) has a similar 
asymptotic expansion 



1 



1 (2.) 

n^ c (n,r 12 ) = n 2 , c (ri,ri 2 ) + -^n 2 '(ri, r X2 ) H . 

p 2 

(C26) 

The absence of term in 1/// in the last expansion im- 
plies that the first two terms in the asymptotic expan- 
sion of dE^/dfi can be determined by considering only 
ri2, c (jc 1, r 12). We find 



dE£ 



VttAo 
M 4 



J ?i 2 , c (r, r)dr 

«2. C ( r : r ) dr + 



(C27) 



where n 2 c (r, r) is the spherical-averaged derivative of the 
on-top correlation pair density which, according to the 
electron-electron cusp condition [6(j , is equal to the total 
on-top pair density: n' 2c (r, r) = n 2 (r,r). Therefore, for 
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large \x the short-range correlation energy has the exact 
behavior 



m = ^J n 2 ,c(r, r)rfr + h&*L J n 2 ( r , r )dr + ■■■. 

(C28) 
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